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Experimental measurements of pressure distributions and force coefficients 
obtained from a squeeze film damper test rig executing a circular centered orbit 
are presented. The test rig has been designed to study the effect of fluid inertia 
on the pressure field and dynamic force response on a damper configuration 
with a relatively large clearance. Past measurements of the squeeze film damper 
force characteristics have been carried out at squeeze film Reynolds numbers 
not exceeding a value equal to 10. In the present paper, following contemporary 
applications, operations at Reynolds numbers up to fifty are tested for CCO’s 
with an orbit radius=0.8 (Re < 10 at e = 0.5). The results obtained from a 
fully open ended damper are presented in detail. The effects of fluid inertia, 
cavitation and the open end geometry on the pressure distributions and force 
coefficients are discussed. 
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NOMENCLATURE 

=SFD radial clearance 

=pu>LR 3 /C damping coefficient conversion factor 
=pLR 3 /C inertia coefficient conversion factor 
—pujR 3 L/C 2 force conversion factor 
=puiR 2 /C 2 pressure conversion factor 
=C P ■ Re = pu> 2 R 2 pressure conversion coefficient 
=dimensionless direct damping coefficients 
normalized by p,R 3 L/C 3 
=damper journal diameter 
=dimensionless direct inertia coefficients 
=f r /e, normalized by pR 3 L/loC 3 

=dimensionless radial and tang, film forces normalized by Cf 
=Jpcos0dO, f p sin QdQ, respectively 
=damper journal length 
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=dimensionless pressure normalized by C p 
=damper journal radius 
=tvC 2 /v modified Reynolds number 
=axial coordinate 
=dimensionless orbit radius 
=circumferential coordinate 
=absolute viscosity 
=kinematic viscosity 
=density 

=90° + tan -1 f r / ft, phase angle 
=frequency of damper journal center motion 

INTRODUCTION 

Squeeze film dampers have been successfully used to control vibration in 
turbomachinery. Numerous studies both theoretically and experimentally have 
been presented for squeeze film dampers. However, due to the fluid film pressure 
distortion factors such as oil grooves, end seals, cavitation, and fluid inertia, 
experimental results are often not consistent with predictions. 

In the 1970’s, many researchers [1-5] measured the hydrodynamic pressure 
field and damping; coefficients from a number of squeeze, filmdamper. test rigs and 
compared them to theiconventionalisolution of Reynolds equation. Correlation be- 
tween experimental measurements and theoretical predictions was qualitatively 
reasonable but quantitatively poor in most cases. 

In the 1980’s, the effect of fluid inertia on the pressure distribution and 
dynamic force response of squeeze film dampers was shown to be of great impor- 
tance [6-10]. Experimental measurements [11-15] demonstrated the significant 
effect of fluid inertia through both the direct measurements of film pressure 
and forces and the indirect measurements from the dynamic response of rotor 
systems. 

Most measurements of a squeeze film damper have been carried out at 
squeeze film Reynolds numbers not exceeding 10. However, in contemporary 
practical applications, the squeeze Reynolds number ranges from 10 to 50. Thus 
the need exists for experimental measurements to cover the maximum practical 
range of Reynolds numbers up to Re=50. 

The objective of this paper is to present an experimental study on the effect 
of fluid inertia and cavitation on the pressure distributions and force coefficients. 
Two types of squeeze film dampers with its journal executing a circular centered 
orbit are tested: fully open ended and partially sealed configurations. The results 
obtained from a fully open ended test rig and a partially sealed test rig are 
presented in part I and II, respectively. The range of Reynolds number tested 
is between 2 and 50. Measured fluid film pressure profiles and force coefficients 
axe compared with existing theoretical predictions. 
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EXPERIMENTAL FACILITY AND PROCEDURE 

Description of the SFD Test Apparatus 

The schematic drawing in Figure 1 shows the structure of the squeeze film 
damper (SFD) test rig. The damper journal with an outer diameter of 12.7 
cm (5 inch) and a length of 2.4 cm (0.94 inch) is mounted with a press fit on 
the outer race of a ball bearing. An O-ring seal is installed at the left end 
face of the damper journal and provides a full proof of a no leak seal. The L/D 
ratio is equal to 0.188 and simulates a short bearing configuration typical of a 
real engine damper application. The clearance is 0.159 cm (0.0625 inch). The 
large clearance produces large values of Reynolds number at low speeds so that 
significant effects of fluid inertia axe produced. 

The end plate of the journal secures the right end face of the journal. The 
outer diameter of the end plate is used to measure the orbital motion of the 
journal during operation. The end plate not only preserves alignment of the 
housing and journal, but also provides a circumferential end cavity. The width 
and depth of the end cavity axe 0.44 cm (11/64 inch) and 0.40 cm (0.156 inch), 
respectively. The end cavity is utilized with ancillary hardware to produce 
several different configurations such as short, long, and finite length bearings 
using different types of end seals. 

The ball bearing supporting the damper journal is mounted on a replaceable 
eccentric sleeve that lockes to the shaft with a special retainer. The shaft is 
driven by a constant speed electric motor at 1770 rpm (29.5 Hz). The motor 
end of the shaft is supported by two more ball bearings (not shown in Figure 
1) mounted on stiff supports to minimize lateral movement of the shaft and the 
damper journal. The shaft is connected to the motor by a flexible coupling to 
prevent axial loads on the journal. 

The bearing housing has numerous holes for oil inlet and outlet, pressure 
transducers, and proximity probes. Figure 2 shows a sketch of the squeeze film 
land. At the axial location Z0, two 0.24 cm (3/32 inch) diameter holes located at 
the circumferential locations 270° (top) and 90° (bottom) serve as the oil inlets 
to the damper annulus. The two oil inlet holes have check valves to prevent flow 
reversals from the squeeze film region. Dial pressure gauges are connected to 
the inlet lines in order to measure the oil supply pressure. 

At two different axial locations, Zl=0.56 cm (0.22 inch) and Z2=1.67 cm 
(0.66 inch), a total of 16 holes axe distributed in the circumferential direction for 
installation of pressure transducers (PT). These holes may be used to directly 
measure the temperature of the fluid film inside the damper using a dummy 
pressure transducer with a small hole plugged with a thermocouple wire. The 
unused holes in the housing axe plugged with dummy pressure transducers to 
prevent lubricant leakage. 

At the axial location Z3=3.72 cm (1.46 inch), a single 0.64 cm (1/4 inch) 
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diameter oil outlet is located at the circumferential location 270° (top) of the 
housing in order to get a completely submerged arrangement. Oil is never 
pumped out of the damper oil outlet, but flows naturally. The bearing housing 
has a hole for proximity probes (PP) at the axial location Z3 and circumferential 
locations 180° (left). 

The Plexiglas cover with a lip seal covers the end of the bearing housing 
and insures a completely submerged arrangement and sufficient oil in the squeeze 
film region. The completely submerged arrangement prevents air entrance into 
the squeeze film region. 

Experimental Procedure 

Figure 3 shows the instrumentation arrangement for measurement of the 
fluid film pressure and real time calculation of the film force coefficients at a 
specified axial location of the SFD test rig. 

In the case of a circular centered SFD, since the time variations in cir- 
cumferential pressure measured at any location on the damper annulus should 
be identical for any one cycle of the journal, the fluid film pressure profile for 
a centered orbit motion can be measured with only one pressure transducer. 
Furthermore, since the orbit motion is centered, only one proximity probe is 
required to determine the dynamic eccentricity ratio. 

Therefore, one proximity probe was installed at the axial location Z3 and 
the circumferential location 180° on the outside diameter of the end plate. The 
proximity probe measures the journal displacement, from which the variation 
of film gap in the squeeze film region can be determined. Two dynamic 
pressure transducers were installed at two different axial locations Z1 and Z2 
and circumferential location 180° to study the axial pressure variation caused 
by end leakage. The variation of the pressure distribution in both the axial 
and the circumferential directions are monitored with two dynamic pressure 
transducers. The piezo-electric pressure transducers used can measure only the 
dynamic variation of film pressure. Thus, a strain-gage pressure transducer 
to measure absolute pressure was installed at the axial location Z1 and the 
circumferential location 330°. Absolute value of pressure obtained from the 
pressure transducer gives information about the cavitation pressure. 

A thermocouple type T was embedded with an epoxy resin into a small 
hole drilled into a dummy pressure transducer at the axial location Z1 and the 
circumferential location 210°. The thermocouple is flush with the wall of the 
bearing housing so that the temperature of the oil film in the squeeze film region 
is directly measured. The value of temperature was input for calculation of the 
viscosity of the oil and the corresponding Reynolds number. The increase of 
the temperature of the oil is mainly due to shear friction of the oil and heat 
conduction from the ball bearing supporting the damper journal. 

The squeeze film Reynolds number Re — loC 2 jv physically represents the 
ratio of the fluid inertia force to the viscous force. Observation of the variation 


82 



in Reynolds number is very important to study the effects of fluid inertia on 
the pressure distribution and force coefficients. Since the whirling frequency u> 
of the journal and the radial clearance C of the SFD test rig are fixed values, 
the oil viscosity is the only factor that can be varied to obtain a change in the 
Reynolds number. 

Oil viscosity was measured fora rangeofitemperatures prior to experimental 
tests using a Synchro-Lectric viscometer (LVF model). Oil viscosity formulae 
were obtained in terms of temperature using a mathematical relation suggested 
in the ASTM D-341. Three kinds of oils have been used in the experimental tests; 
SAE30 oil (engine oil), IS022 oil (turbine oil), and DOWIO oil (silicone oil). 
The range of temperature measured was between 27° (7(80° F) and 49° <7(120° F). 
The corresponding Reynolds numbers range from Re=2 to Re=7 for SAE30 oil, 
Re=17 to Re=35 for IS022 oil, and Re =42 to Re=48 for DOWIO oil. 

Oil supply pressures were about 68.75 KPa (10 Psig) which allowed a flow 
rate equal to 200 ml per min. for SAE30 oil, and about 13.8-20.6 KPa (2-3 Psig) 
for lower viscosity oils. Higher oil supply pressures can not be applied to the 
damper because of deformation of the plexiglas cover. Also, higher oil supply 
pressures produce a notable distortion of the fluid film pressure field, because 
the two small oil inlets directly supply oil into the squeeze film region (there 
is no circumferential oil inlet groove). The oil inlet groove is often employed 
in squeeze film dampers to prevent distortion of the pressure field due to inlet 
holes. Previous experimental measurements have shown that an oil groove has a 
large effect on the pressure distortion [11,13]. Since experiments with a circular 
centered orbit should yield pressure distributions independent of the location of 
the pressure transducer, it is very important to get rid of potential distortion 
factors for better comparison between experiment and analysis. 

In the experimental measurement reported here, 60 data points per cycle 
per pressure transducer and proximity probe were taken. Each measurement 
consisted typically of 3 cycles. The force coefficients are the averaged values 
obtained from integration of these 3 cycles of the measured pressure profiles. 

EXPERIMENTAL RESULTS AND DISCUSSIONS 

Experimental tests have been performed on the squeeze film damper (SFD) 
test rig to measure dynamic pressure distributions and force coefficients with 
orbit radii equal to 0.5 and 0.8. Table I shows the conversion factors for proper 
dimensional values of pressure and force coefficients for the three kinds of oils 
used in the experimental work. These factors are used to get the corresponding 
dimensional values of the measured dimensionless pressure and force coefficients. 
To predict the measured pressure distributions and film force coefficients, a 
modified short bearing solution [16] is used. The analytical model used included 
the effects of finite length and fluid inertia. 
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Pressure Distributions 

Figures 4 and 5 show comparisons between the measured and predicted 
pressure profiles for orbit radius equal to 0.5 and 0.8, respectively. The viscous 
lubrication solution produces antisymmetric pressure profiles with respect to the 
maximum film thickness (u;2* = 7r), while the additional pressure due to the fluid 
inertia is symmetric with respect to cut* = i r [14]. Referring to the combined 
pressure profile including the additional inertial pressure, the main effects of fluid 
inertia are that the positive dynamic peak pressure is smaller than the negative 
dynamic peak pressure, and the dynamic pressure has a negative value at the 
minimum film thickness (cut* = 0) while it has a positive value at the maximum 
film thickness (cut* = 7 r). 

The pressure profile with a orbit radius e = 0.8 corresponds to a cavitated 
film as shown in Figure 5. The absolute levels of pressure measured are close to 
zero absolute and thus, these indicate vapor cavitation in the lubricant film. 

According to conventional lubrication theory, negative dynamic pressure is 
produced in the diverging section of bearing while positive dynamic pressure is 
produced in the converging section of bearing. Contrarily, as shown in Figure 
5, negative pressure exists in the converging section of the SFD. It has been 
generally assumed that the positive pressure profile is unaffected by cavitation, 
but it is presumed here that vapor cavitation influences the positive pressure 
region. 

For reduced oil viscosities, the extent of the vapor cavitation becomes 
smaller. For e = 0.8, the measurements showed that the cavitation completely 
disappeared at Re=22.1. Figure 6 shows the measured pressure which corre- 
sponded to IS022 oil. 

An important characteristic of the open ended configuration is the large 
pressure drop in the axial direction due to the substantial flow in the axial 
direction. In Figures 5 and 6, the amount of pressure drop measured at axial 
location Z2 is about 40% of the peak positive pressure at axial location Zl. 
The modified short bearing solution [16] predicts 46% of pressure drop between 
the two axial locations. The measurements confirm that the modified short 
bearing model can adequately describe the axial pressure drop in a parabolic 
form. However, the predicted pressure profile at the axial location Z2 is poor 
in the negative dynamic pressure region. The peak negative pressure drop 
measured at the axial location Z2 is about 20% of that at Zl at large Reynolds 
numbers (Re>20). The parabolic variation of the axial pressure predicted 
from the modified short bearing solution is therefore not valid in the negative 
pressure region where the Reynolds number is large. Since the film is completely 
submerged, the negative pressure produced in the diverging section creates 
reversed flow in the axial direction. Thus, a Bernoulli effect is likely to be 
significant which violates the boundary condition p— 0 at the end of the damper. 

Figure 7 shows dynamic pressure waves obtained at Reynolds number equal 
to 44.2. In these figures the distortion of the antisymmetric shape of the viscous 
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pressure profile is magnified, and the pressure profile is dominated by the effects 
of fluid inertia. The negative peak pressure is about three times the positive 
peak pressure in magnitude. The analytical predictions deteriorate for these 
large Reynolds numbers. 

During the pressure measurements at small Reynolds number with e— 0.5 
and at large Reynolds number with e=.8, the pressure signals on an oscilloscope 
were fluctuating, as seen on the Figures 4 and 7. Thus the fluid flows at these 
large Reynolds numbers may not be laminar and the Bernoulli effects may 
become more significant at low viscosities. Also there is a noise interference 
at low levels of pressure. For an orbit radius e = 0.5, it was not possible to 
measure the pressure distribution at large Reynolds number (Re>10) because 
of the large ratio of noise to measured values. 

During operation of the SFD test rig no bubbles came out of the squeeze 
film region, but many minute bubbles simultaneously appeared in the oil a few 
seconds after stopping the test rig and oil supply pump. It is presumed that 
the air dissolved in the oil came out of solution when the external oil supply 
pressure was removed. 

Force Coefficients 

Experimental fluid film forces and force coefficients are determined by 
numerical integration of the measured pressure waves around the damper journal 
surface. For the analytical predictions, a full film assumption (2tt film) with 
orbit radii e = 0.5 and 0.82 are used. Since the pressure profiles for e = 0.8 are 
cavitated at small Reynolds numbers (Re< 10), corresponding predictions axe 
excluded here. The force coefficients are calculated locally and corresponded to 
a unit axial width of the damper at the specific axial locations. 

Figures 8 show the local damping coefficients C u computed from the 
measured pressure profiles. Since cavitation is present in the pressure profiles 
of e = 0.8 and for Reynolds numbers less than 10, the measured damping 
coefficients must be regarded as equivalent damping coefficients, i.e. equal to 
(Cu — D tr )• But at Reynolds numbers larger than 15, the damping coefficients 
axe purely direct coefficients because of no cavitation in the measured pressure 
profiles. 

Generally, the predictions for the local damping coefficients using the 
modified short bearing solution [16] compare favorably with measurements at 
the axial location Zl. However, the predicted coefficients axe poor at the axial 
location Z2. The main reason is that the parabolic axial pressure profile obtained 
from the short bearing theory does not properly predict the axial pressure drop 
in the negative region when both Reynolds number and orbit radius are large. 
A level of the dimensional values of damping achieved in the test rig can be easily 
obtained if the measured damping coefficients Ctt are multiplied by CdRe / Re 
( CdRe is given Table I). 

Figures 9 show the local inertia coefficients D rr calculated from the mea- 
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sured pressure data. For e = 0.8 at Reynolds number below Re=10, the inertia 
coefficients are to be regarded as an equivalent inertia coefficient, i.e. equal to 
(D rr — C rt ) due to the observed cavitation in the measured pressure profiles. 
Since the region of cavitation increases as the Reynolds number decreases, the 
equivalent inertia coefficients correspond mainly to the cross coupled damping 
coefficients at small Reynolds number with e=0.8. The dimensional values of 
inertia coefficients can be obtained by multiplication of /Re in Table I. 

For a orbit radius e = 0.8 at Reynolds number less than Re=10, the inertia 
coefficients at both the axial location Z1 and Z2 are negative, showing the cross 
coupled damping force is dominant. (Cross coupled damping has been described 
as dynamic stiffness [17,18]). For Reynolds number larger than 15, the inertia 
coefficients become positive, which means that fluid inertia overcomes the effect 
of the cross coupled damping and acts as an added mass on the journal. The 
prediction of the local inertia coefficients are poor at the axial location Z2 near 
the end of the damper journal as in the case of the damping coefficient prediction. 
Once again, the parabolic axial pressure profile can not be used in the negative 
pressure region. 

Damping coefficients are seen to be large at small Reynolds number and 
indicate a dominance of viscous effects. As the Reynolds number increases, the 
inertia coefficients eventually approach the same level of the damping coefficients 
as shown in Figures 8 and 9, and show that the fluid inertia forces become 
dominant 

San Andres [9,10] determined fluid film forces and force coefficients based 
on the assumption that fluid inertia do not affect the velocity profiles or flow 
rates. The important results were 

1) The dimensionless damping coefficient Ctt is purely viscous and not 
affected greatly by fluid inertia. 

2) The fluid inertia force increases linearly with Reynolds number so that 
the dimensionless inertia coefficient D rr /Re has a constant value for small to 
moderate Reynolds number. 

However, since the above assumption is valid for very small Reynolds 
number (Re<l), as the fluid inertia forces approach the same order of magnitude 
as viscous forces for large Reynolds number (Re>10), the effect of fluid inertia 
forces on the velocity profile then becomes significant and needs to be considered. 
When the effect of the fluid inertia forces on the velocity profiles or flow rates is 
included in the analysis, the dimensionless damping coefficient C t t and inertia 
coefficient D rr /Re, both become strongly dependent on the Reynolds number 
[16]. Experimental evidence supporting these results is shown in Figures 8 and 
9 . 

Figure 10 shows the force phase angle calculated from the force coefficients 
given in Figures 8 and 9. It is important to notice how the point of application of 
the total fluid film force relative to the position of the maximum film thickness 
changes as the Reynolds number increases. If the viscous film force is dominant, 
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the radial film force is very small compared to the tangential film force and the 
phase angle approaches 0 = 90°; while if the inertial film force is dominant, the 
radial film force is very large (outward) compared to the tangential film force 
and the phase angle <j> approaches <f> = 180°. Thus, phase angles above 90° 
indicate a film radial force component due to fluid inertia. Predictions of phase 
angle obtained from the modified short bearing solution agree well with the 
measurements. 


CONCLUSIONS 

Experimental measurements have been performed on an open ended squeeze 
film damper executing circular centered orbits. Dynamic pressure measurements 
were obtained for two orbit radii, e = 0.5 and 0.8. Corresponding force coeffi- 
cients were calculated by numerical integration of the measured pressure profiles 
From the experimental measurements presented, the following conclusions can 
be described: 

1) The parabolic axial pressure profile obtained from the modified short 
bearing solution can adequately predict the axial pressure drop in the positive 
dynamic pressure region, but not in the negative dynamic pressure region. These 
effects are more pronounced at large Reynolds numbers. 

2) Cavitation observed in the submerged damper corresponds to vapor 
cavitation. The vapor cavitation influences not only the negative pressure field 
in the diverging region of a damper, but also the positive pressure field in the 
converging region of the damper. 

3) Inertia coefficients D rr /Re (radial film force) become comparable to 
damping coefficients C tt (tangential film force) as Reynolds number increases. 
Both the measured inertia and damping coefficients increase with increasing 
Reynolds numbers. 

The experimental results obtained from a partially sealed configuration are 
presented in Part II. 
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conversion t 
factor 
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CpR* 

Kpa 

(P«) 

N sec/m 
(lbf sec jin) 

Cifu 

N sec 2 /m 
(lbf jcc 2 fin) 

SAE30 ' 

1 

119.75 

616.99 

3.328 

1 

(17.37) 

(3.523) 

(0.0190) 

IS022 

i 

123.06 

634.03 

3.420 

i 

(17.85) 

(3.620) 

(0.0195) 

DOWIO i 

I 

128.98 

664.56 

3.584 

t 

1 

(18.71) 

(3.795) 

(0.021) 



TABLE I. Conversion factors for a SFD test rig at 26.7° C 


oil 

outlet 

r 


1 . Housing 

2. End cavity 

3. Journal 

4. Eccentric sleeve 

5. Sleeve holder 

6. End plate 

7. Plexiglas cover 
and Lip seal 

8. Housing support 

9. Shaft 

10. Anti rotation ear 

11. O-ring 


PT: Pressure Transducer 
PP: Proximity Probe 
PQ: Pressure Gauge 
CV: Check Valve 


Fig. 1 Schematic drawing of a SFD test rig 
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SFD geometry of an open ended configuration Fig. 3 Schematic diagram of instruments 
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Dimensionless pressure profiles, Fig. 5 Dimensionless pressure profiles, 

open ended configuration, Re=2.14, e = 0.5 open ended configuration, Re=2.26, 




93 


Dimensionless damping coefficients vs. Reynolds no. pig. 9 Dimensionless inertia coefficients vs. Reynolds no. 
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